Anomalous roughening of Hele-Shaw flows with quenched disorder 
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The kinetic roughening of a stable oil-air interface, moving in a Hele-Shaw cell which contains 
a quenched columnar disorder (tracks) has been studied. A capillary effect is responsible for the 
dynamic evolution of the resulting rough interface, which exhibits anomalous scaling. The three 
independent exponents needed to characterize the anomalous scaling are determined experimentally. 
The anomalous scaling is explained in terms of the initial acceleration and subsequent deceleration 
of the interface tips in the tracks coupled by mass conservation. A phenomenological model that 
reproduces the measured global and local exponents has been introduced. 
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The kinetic roughening of surfaces and interfaces in 
nonequilibrium conditions has received a great deal of 
attention in last years. This is a subject of technolog- 
ical importance in surface growth, and of fundamental 
relevance as an example of scaling far from equilibrium 

Roughening originates from the competition between 
stabilizing and destabilizing mechanisms at different 
length scales. As a result of the competition, the cor- 
relation length of the interfacial fluctuations along an 
initially flat front increases monotonously with time, un- 
til it reaches the system size and the roughness satu- 
rates. The statistical properties of the roughening front 
are usually given in terms of the (rms) width w of the 
interfacial fluctuations. In the absence of characteristic 
scales, many systems follow the dynamic scaling hypoth- 
esis of Family-Vicsek (FV) §, in which short and long 
length scales have the same scaling behavior: w(l) ~ 
for t < t x , ~ l a for t > t x and t x ~ I*. Here I is the 
lateral size of the window used to measure w, t x is the 
saturation time, (3 the growth exponent, a the roughness 
exponent, and z the dynamic exponent. The scaling re- 
lation a — z[3 reduces to two the number of independent 
exponents that characterize the universality class of the 
problem under study. 

In recent years, however, experiments on molecular 

beam epitaxy ||, sputtering [Ej, fracture mechanics ||, 

and electrodeposition || have revealed a rich variety of 

situations in which the scaling of the global and the local 

interface fluctuations are substantially different. Long 

length scales are then characterized by the same FV 

scaling exponents, while short length scales follow the 

so-called anomalous scaling ansatz [Q: w(l,t) ~ t 13 for 

t < t t , w(l,t) ~ l a i°°tP* for t t <t <t x , and w{l,t) ~ 
ia lBcLa - aioe for tx < t The local time tl ^ and the 



saturation time t x ~ L z , where L is the system size in the 
direction along the front. Here (3* and a/ oc are new lo- 
cal growth and roughness exponents, respectively, which 
characterize the scaling at short length scales (I <C L). 
There is a new scaling relation a — ai oc — z(3* . Since 
now there are five exponents and two scaling relations, 
the anomalous scaling is characterized by three indepen- 
dent exponents. Usually the roughness exponents are 
determined through the power spectrum, which is less 
sensitive than w to finite size effects. The spectrum 
scales as S(k,t) = /c^( 2a + 1 ) s^fci 1 / 2 ), where sa obeys 
sa{u) ~ u 2a+1 when u 1 and sa(u) ~ u 29 when 
U ^> 1. Here = a — ai oc . The anomalous scaling 
leads to the usual FV scaling when = (short and 
long length scales are characterized by the same rough- 
ening exponents) and 0* = (the temporal correlations 
of the short length scales disappear). Recently it has been 
shown that the anomalous scaling can be related to the 
presence of high slopes in the front that have nontrivial 
dynamics [p|. 
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FIG. 1. Sequence of typical oil-air interfaces. The sili- 
cone oil moves upwards in the picture, and the disorder pat- 
tern is represented in grey. The experimental parameters are 
b = 0.36 mm and v = 0.04 mm/s. 
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In this Letter we report the first observations of anoma- 
lous roughening in experiments of Hele-Shaw flows with 
quenched disorder. We consider an initially flat front 
between oil and air, driven at a small constant average 
velocity v in a horizontal Hele-Shaw cell with a columnar 
quenched disorder. The disorder consists of continuous 
copper tracks on a fiberglass substrate in the advanc- 
ing direction of growth and randomly distributed in the 
perpendicular direction (Fig. |l|). In this situation, the 
correlation of the disorder in the advancing direction is 
infinite. The local motion relative to the average position 
of the front is driven by capillary forces, caused by the 
different curvatures of the advancing front in the third 
dimension, depending on whether the oil is on a copper 
track or on the fiberglass substrate. An analysis of the 
rms width of the rough interface at short and long length 
scales plus a collapse of the power spectra of the interfa- 
cial fluctuations enables an independent determination of 
the three exponents characteristic of anomalous scaling. 

We associate the anomalous roughening with the initial 
acceleration and subsequent deceleration of the liquid on 
the copper tracks, caused by capillarity, and the coupling 
of the motion over tracks and over fiberglass due to mass 
conservation. The anomaly decreases as v increases, and 
disappears when v is large enough to override the local 
acceleration and deceleration To explain the exper- 
imental results, we propose a phcnomenological model 
which gives global and local exponents in good agree- 
ment with the experimental exponents. 

The experiments have been performed in a horizon- 
tal Hele-Shaw cell, 190 x 550 (L x H) mm 2 , made of 
two glass plates 20 mm thick. A fiberglass substrate 
with copper tracks, manufactured using printed circuit 
technology, is attached to the bottom plate. The unit 
copper track is d = 0.06 ± 0.01 mm high and 1.5 ± 0.04 
mm wide. The tracks are placed randomly in a grid of 
1.5 mm lattice spacing, along the x direction, up to a 
35% occupation. There are no lateral correlations of the 
disorder at distances larger than the width of the unit 
track. The distance between the top plate and the sub- 
strate defines the gap spacing b, which has been varied 
in the range 0.16 < b < 0.75 mm. We have used 4 dif- 
ferent disorder configurations and carried out 2 identical 
runs per disorder configuration. A silicone oil (kinematic 
viscosity v — 50 mm 2 /s, density p — 998 kg/m 3 , and 
surface tension oil-air a = 20.7 mN/m at room temper- 
ature) is injected into one side of the cell at constant 
volumetric injection rate. The oil wets similarly and al- 
most completely the surfaces which it is in contact with 
(glass, substrate, and copper). The evolution of the in- 
terface is monitored using two CCD cameras mounted 
along the growth direction and fixed. The final resolu- 
tion is 0.37 x 0.37 mm 2 /pixel, providing 4 pixels per unit 
track in the x direction. We have carried out a series 
of reference experiments to calibrate side walls effects 
and possible inhomogeneities in gap spacing. Our results 



make us confident that these disturbances do not have a 
relevant effect on the dynamic scaling presented in this 
work. Extensive details of the experimental setup will be 
given in a forthcoming publication |ic|] . 
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FIG. 2. Experimental determination of (3 and /3* . The win- 
dow size I decreases from L to L/256. The experimental pa- 
rameters are b = 0.36 mm and v = 0.08 mm/s. 



In Fig. |2j we present the results for the interfacial width 
as a function of time at different length scales /, from L to 

L/256. We plot W(l,t) = [w 2 (l, t) - w 2 (l, 0)]^ 2 to min- 
imize the influence of the initial condition |flpT| . We ob- 
tain clear power law dependencies, providing (3 = 0.52 ± 
0.02 at the longest length scales and (3* = 0.26 ± 0.03 
at the two shortest length scales. It is worth mentioning 
that a different method, based on the growth of the local 
slopes § , gives the same value of (3* . The fact that this 
local exponent is significantly different from is a sign of 
anomalous scaling of the measured fronts. 

The temporal evolution of the power spectrum is pre- 
sented in Fig. ^(a). At very short times we observe a 
regime with S(k) ~ fc~ 3 ' 4 . It is a super-rough tran- 
sient regime created when the fluid on the tracks reaches 
its maximum velocity after coming into contact with the 
disorder. The main regime appears further on, with 
S(k) ~ k~ , and a vertical shift that progressively de- 
creases with time and disappears at saturation (t > 450 
s) . The vertical shift of the power spectra is a sign of in- 
trinsic anomalous scaling, in which the spectral exponent 
JjJI can be identified with the local roughness exponent, 
so that ai oc = 0.55 ± 0.10. 

Using the above procedure, the three independent ex- 
ponents (3, (3* , ai oc have been obtained directly from the 
experiments. The other two exponents a, z can be ob- 
tained from the scaling relations. In order to verify the 
whole scaling, however, we prefer to determine the crit- 
ical exponents by performing the best collapse of the 
spectra compatible with the experimental results. The 



2 



collapse, presented in Fig. |](b), leads to the following set 
of scaling exponents: 

(3 = 0.50 ± 0.04, 0* = 0.25 ± 0.03, 

a = 1.0 ± 0.1, a loc = 0.5 ± 0.1, z = 2.0 ± 0.2. (1) 
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order q, of the form {(h(x + l,t) - h(x, t)) 9 } 1 / 9 - l a « . 
Different exponents for a q : a.2 — 0.6, ~ 0.3, and 
ag ~ 0.2 were obtained, which confirm the existence of 
multiscaling in the experiment. 

To understand the origin of anomalous scaling, we 
should look at the interfacial dynamics in greater de- 
tail. An initially flat oil-air interface, moving at a nom- 
inal average velocity v, experiences local accelerations 
at those points that touch a copper track for the first 
time. The local velocity of these points jumps to a max- 
imum in a characteristic time of 1 s, and then relaxes as 
t~ x / 2 to v. The maximum velocity reached in a copper 
track increases with track width and decreases with gap 
spacing. As shown in Fig. |], the average interface veloc- 
ity over copper tracks or over fiberglass channels follows 
v± = v ± (vm — w)t -1 / 2 respectively, where vm is the 
maximum of the average velocity over tracks. This func- 
tionality varies slightly with gap spacing due to the in- 
creasing correlation between neighboring tracks. It does 
not apply for b > 0.6 mm. This behavior of different 
velocities v+ and u_ is directly correlated with the shift 
of the power spectra, indicative of anomalous scaling: 
the shift is maximum at very short times and disappears 
close to saturation, similarly to the difference between 
the average velocities shown in Fig. |l|. When we inhibit 
the relaxation of the velocity over copper tracks by ei- 
ther injecting at v > Vm or using large 6, the anoma- 
lous scaling disappears. This scenario has been studied 
experimentally by exploring different velocities and gap 
spacings. Fig. ^| shows a phase diagram representing the 
regions where the anomalous scaling is present (grey re- 
gion, 8 ^ 0) or not (white regions, 9 — 0). The solid line 
represents the function and the vertical line the 

limit between strong and weak capillary forces. 
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FIG. 3. (a) Temporal evolution of the power spectrum for 
the same experiments as Fig. ^. The vertical line gives the 
value of k associated with the lateral size of the unit cop- 
per track, (b) Collapse of the experimental power spectra for 
t > 10 s. (c) Collapse of the power spectra obtained from 
numerical simulations of the interfacial equation (|^). 

We have also looked into the possibility of multiscaling 
p3| through the scaling of the generalized correlations of 
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FIG. 4. Average interface velocity over copper tracks and 
over fiberglass, for b = 0.36 mm and v — 0.08 mm/s. 



The critical exponents of the anomalous scaling de- 
pend on gap spacing and velocity. We have seen that 
the anomalous exponent 9 tends to 9 = as we increase 
either the velocity or the gap spacing. This exponent 
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is particularly sensitive to variations of the gap spacing. 
For b > 0.6 mm, 8 = at any velocity. The experimental 
parameters b = 0.36 mm and v = 0.08 mm/s, deep in the 
grey region of Fig. ^[ give the appropriate conditions to 
fully characterize the anomalous scaling experimentally. 
Smaller gap spacings mean that the correlation length 
in the x direction cannot grow to scales larger than the 
average track width, and larger gap spacings give too 
weak an anomalous scaling because the t -1 / 2 law is only 
obeyed on the widest tracks. 
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FIG. 5. Phase diagram v(b) indicating the regions where 
the anomalous scaling is observable. The symbols represent 
the different regions explored experimentally, and the arrow 
indicates the parameters used in the paper. 



shown in Fig. ||(c), and the multiscaling exponents, repro- 
duce the values determined experimentally. In summary, 
Eq. (^) enables the relevant effects in the experiment to 
be both identified and calibrated. 

In conclusion, we have presented the first experimental 
evidence of anomalous scaling in the dynamic roughening 
of a fluid interface in a Hele-Shaw cell with quenched dis- 
order, and we have studied the experimental conditions 
for the appearance of the observed anomaly. The expo- 
nents (3, (3* , and a\ oc have been measured independently 
and their value has been refined by imposing the scaling 
relations through a collapse of the power spectra. Finally 
we have introduced an interfacial equation (|J) that mod- 
els the capillary phenomena observed in the experiments, 
at a phenomenological level. The model reproduces both 
the morphology of the interfaces and the values of the 
anomalous scaling exponents. 
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Finally, in order to reproduce the experimental behav- 
ior in a phenomenological way, we have developed an 
interfacial equation of the form: 

d t h = d x D{x)d x h + v + (v M - v)t~ 1/2 v(x), (2) 

where h = h(x,t), and rj(x) is a dichotomous noise with 
values +1 and —1 (columnar disorder). The diffusive 
coupling accounts for the dynamics with z — 2 found in 
experiments, and the other two terms account for the be- 
havior of v± observed experimentally. Although the term 
t~ V 2 introduces some nonlocality in the model, this is es- 
sentially different from non-local models |L4j introduced 
recently to account for mass conservation. The simplest 
version (constant D) reproduces the experimental behav- 
ior at long length scales (a = 1, (3 = 0.5, z = 2), but not 
at short length scales, since numerical simulations of the 
model give (3* — and ai oc = a = 1. The reason is that 
the simulated interfaces do not display the large fluctu- 
ations at the copper track edges observed in the experi- 
mental fronts (Fig. |l|). This can be imposed analitically 
in Eq. (|^) through an inhomogeneous diffusion coefficient 
D(x). Numerically, we have introduced this effect by spa- 
tially averaging the interface inside each track once each 
n time steps. We thus recover interfaces that are morpho- 
logically analogous to experimental interfaces. Now the 
exponents derived from the scaling of the power spectra, 
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